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Defect generation and deconfinement on corrugated topographies 
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We investigate topography-driven generation of defects in liquid crystals films coating frozen 
surfaces of spatially varying Gaussian curvature whose topology does not automatically require 
defects in the ground state. We study in particular disclination-unbinding transitions with increasing 
aspect ratio for a surface shaped as a Gaussian bump with an hexatic phase draped over it. The 
instability of a smooth ground state texture to the generation of a single defect is also discussed. 
Free boundary conditions for a single bump are considered as well as periodic arrays of bumps. 
Finally, we argue that defects on a bump encircled by an aligning wall undergo sharp deconfinement 
transitions as the aspect ratio of the surface is lowered. 

PACS numbers: Valid PACS appear here 



I. INTRODUCTION 

The melting of a two dimensional crystal can occur 
continuously via two second order topological phase tran- 
sitions characterized by the successive unbinding of dislo- 
cation and disclination pairs. At low temperatures, dis- 
locations are suppressed due to their large energy cost, 
but as the temperature is increased, the entropy gained 
by creating defects overcomes their cost in elastic energy 
and dislocation unbinding occurs to reduces the over- 
all free energy of the system [H Qj Q- The quasi- long 
range order of the crystal is thus destroyed leading to 
an hexatic phase that still preserves quasi-long range ori- 
entational order. This phase can be characterized by 
a complex order parameter with six-fold symmetry. As 
the temperature is increased still further, an additional 
disclination-unbinding transition occurs and the hexatic 
order is finally lost in an isotropic liquid phase jl|. 

Experimental evidence for hexatic order and defect- 
mediated melting has been obtained in systems as di- 
verse as free standing liquid crystal films Q, Langmuir- 
Blodgett surfactant monolayers p, two-dimensional 
magnetic bubble arrays [a, electrons trapped on the sur- 
face of liquid helium 0, la 13 , two-dimensional colloidal 
crystals [iflllll and self- assembled block copolymers d^. 

The unbinding of defects in the plane is entropically 
driven and at low temperature defects are tightly bound. 
By contrast, on surfaces with non zero (integrated) Gaus- 
sian curvature, excess defects must be present even at 
very low temperatures. The theory of topological defects 
in ordered phases confined to frozen topographies with 
positive or negative Gaussian curvature has been investi- 
gated previously; see, e.g., [l^i Ej, Eal • As a general rule, 
regions of positive or negative curvature (valleys, hills 
or saddles) lead to unpaired disclinations in the ground 
state, possibly screened by clouds of dislocations. These 
clouds can in turn condense into grain boundaries at low 
temperature. The predictions of recent studies of crys- 
talline order on a sphere 16] have been confirmed in el- 
egant studies of colloidal particles packed on the surface 
of water droplets in oil [I3| ■ Investigations of the physics 
of defects in curved spaces have also been carried out for 



fluctuating geometries [TM Il9l 12(1 l2ll | . The dynamics of 
hexatic order on fluctuating spherical interfaces was stud- 
ied in Ref. ^22] . Quenched random topographies in the 
limit of small deviations from flatness were investigated 
in Ref. [H . 

In the present work, we investigate topography-driven 
generation of defects on simple frozen surfaces with spa- 
tially varying Gaussian curvature whose topology does 
not automatically enforce their presence in the ground 
state. We study in particular a two-dimensional "bump" 
with a Gaussian shape and dimension large compared to 
the particle spacing. For such a hilly landscape, flat at in- 
flnity, the geometric control parameter is an aspect ratio 
given by the bump height divided by its spatial extent. 
Consider a hexatic phase draped over such a bump. For 
small bumps, the ideal hexatic texture is distorted, but 
there are no defects in the ground state. As the aspect 
ratio is increased, we find that disclination pairs progres- 
sively unbind at T = in a sequence of transitions oc- 
curring at critical values of the aspect ratio. The defects 
subsequently position themselves to partially screen the 
Gaussian curvature. For bumps embedded in surfaces 
of sufflciently small spatial extent, a second instability 
of the smooth ground state needs to be considered. In 
this case, the energy stored in the field can be lowered 
by generating a single positive defect at the center of 
the bump. Novel effects also arise when the hilly sur- 
face is encircled by an aligning circular wall that insures 
a 27r rotation of the orientational order in the ground 
state. In this case, some of the positive defects required 
to match the curvature of the boundary are confined to 
a hemispherical cup centered on the bump, provided the 
aspect ratio a is larger than a critical value ao- When 
a is lowered below ao, the positive defects originally 
"trapped" in the hemispherical cup start undergoing a 
series of sharp " deconfinement transitions" , as they pro- 
gressively migrate to new equilibrium positions dictated 
by boundary conditions and the finite system size. We 
also suggest possible ground states for periodic arrays of 
bumps, like those on the bottom of an egg carton. 

A natural arena to experimentally study the interplay 
between geometry and defects is provided by thin copoly- 



mer films on Si02 patterned substrates |23|- Flat space 
experiments by Segalman et al. have already demon- 
strated that spherical domains in block copolymer films 
form hexatic phases |1^ . 

Our results for hexatics on frozen topographies also ap- 
ply to other XY-like models, as might be appropriate for 
tilted surface-active molecules on curved substrates with 
interactions which favor alignment. The results are rele- 
vant as well to two-fold nematic order on frozen topogra- 
phies. In both cases, we expect qualitatively similar de- 
fect unbinding transitions, although the equivalence be- 
comes more exact in the one Frank constant approxima- 
tion |2J| . Related results have been obtained recently for 
order on a torus J23l • Even though the integrated Gaus- 
sian curvature vanishes, defects appear in the ground 
state in the limit of fat torii, unless the number of de- 
grees of freedom is very large. 

The outline of this paper is as follows. In Section II, the 
relevant mathematical formalism is introduced and used 
to highlight similarities and differences between defects 
on surfaces of varying curvature and electrostatic charges 
in a non-uniform background charge distribution in flat 
space. As an example, we calculate the distorted, but 
defect-free, ground state texture of a hexatic confined to 
a surface shaped as a " Gaussian bump" for aspect ratios 
below the first disclination-unbinding instability. In Sec- 
tion III, we investigate curvature-induced defect forma- 
tion for an isolated bump and a periodic array of bumps. 
In section IV, defect deconfinement is discussed and in 
section V various experimental issues related to our anal- 
ysis are highlighted along with some directions for future 
work. The development of the mathematical formalism 
is largely relegated to Appendices. In Appendix 1X1 the 
Green's function for the covariant Laplacian is derived 
by means of conformal transformations. In Appendix B, 
we introduce a geometric potential whose source is the 
Gaussian curvature. In Appendix|21 we present the gen- 
eral formula for the energy of textures with defects in 
terms of the two functions derived in Appendix ^ and 
IbI We thus explore the existence of position-dependent 
defect self-interactions that arise from the varying Gaus- 
sian curvature. Finally, boundary effects are discussed in 
Appendix ^ 



II. HEXATIC ORDER ON A SURFACE 

A. Electrostatic analogy 

The free energy for hexatic degrees of freedom embed- 
ded in an arbitrary frozen surface can be written as 



F 



Ka 
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I dADan^{u)D'^np{vL) , 



(1) 



where u = {ui,M2} is a set of internal coordinates, n(u) 
is a unit vector in the tangent plane, 2?^ is the covari- 
ant derivative with respect to the metric of the surface 
and dA is the infinitesimal surface area llSl Il9l l2ll |2(1 . 



The generalization to systems with a p-fold symmetry 
is straightforward provided that the one Frank constant 
approximation is used and the consequences of the uni- 
axial coupling neglected 1241 . This choice of free energy 
implies that the minimal energy configuration will be 
given locally by neighboring n(u) vectors which differ 
only by parallel transport. The curvature of the surface 
induces "frustration" in the texture. In fact, by Gauss' 
" Theorema egregium" [23, |23| , tangent vectors parallel 
transported along a closed loop are rotated by an amount 
equal to the Gaussian curvature integrated over the en- 
closed area. On a sphere, for example, the hexatic ground 
state always has twelve excess disclinations as a result of 
this frustration [Ij, |23 . More generally, the sum of the 
topological charges on any closed surface is equal to the 
integrated Gaussian curvature. 

By introducing a local bond- angle field 6'(u), corre- 
sponding to the angle between n(u) and an arbitrary 
local reference frame, we can rewrite the hexatic free en- 
ergy introduced in Eq.(^) as: 



F = -K, 



dAg''''{d^e-A^){d0e-Ap), (2) 



where dA = d^u^, g is the determinant of the metric 
tensor g^fs and A^ is the spin-connection whose curl is 
the Gaussian curvature G{u) |23,|23|. The spin connec- 
tion can be viewed as a " geometric vector potential" . A 
free energy like Eq.(2) also describes the charged Cooper 
pairs implicit in the London theory of a superconduc- 
tor well below Tq- In the superconductor analogy, the 
Gaussian curvature plays the role of a (spatially vary- 
ing) external magnetic field. For the problem considered 
here, however, there are interesting new nonlinear effects 
associated with spatial variations in the metric. 

A detailed analysis of the free energy of Eq.(|2Jl for a 
bumpy surface with free and circular boundary condi- 
tions is presented in Appendices (Q and l^)). Here we 
only sketch the main steps and conclusions. The free en- 
ergy can be readily converted into a Coulomb gas model 
by using the relation 

r^d^ {dp0 - A^) = s(u) - G(u) = 7i(u) , (3) 

where 7"^^ is the covariant antisymmetric tensor, G(u) is 
the Gaussian curvature and s(u) = -^ J2i=i Qi^i^ ~ ^i) 
is the disclination density with Nd defects of charge qi at 
positions u^. The final result is an effective free energy 
whose basic degrees of freedom are the defects themselves 

MM 

F ^ ^ [dA jdA' n(u) r(u, u') n(u') , (4) 

where n(u) is defined in Eq.Q. The Green's function 
r(u, u') is calculated (see Appendix^ by inverting the 
Laplacian defined on the surface 



(5) 



r(u,u')^- -T 



and we have suppressed for now defect core energy con- 
tributions which reflect the physics at microscopic length 
scales. Eq.Q can be understood by analogy to two 
dimensional electrostatics, with the Gaussian curvature 
G(u) (with sign reversed) playing the role of a non- 
uniform background charge distribution and the topolog- 
ical defects appearing as point-like sources with electro- 
static charges equal to their topological charge qi. As a 
result, the defects tend to position themselves so that 
the Gaussian curvature is screened: the positive ones 
on peaks and valleys and the negative ones on the sad- 
dles of the surface. However, this analogy does neglect 
position-dependent self- interactions |33, but since these 
are quadratic in the charge they are negligible for hexat- 
ics. Hence positive disclinations of minimal topological 
charge qi = ^ continue to be attracted to positive cur- 
vature (see Appendix 0. 

More generally, we can consider p-fold symmetric or- 
der parameters with minimum charge defects ± — . The 
case p = 1 corresponds to tilt order of absorbed molecules 
and p = 2 describes 2D nematics. The cases p — A and 

M- 6 describe tetradic and hexatic phases respectively 
. Strictly speaking, Eq.(l) only describes the cases 
p = 1 and p = 2 in the one-Frank-constant approxima- 
tion J2J| . Most of our discussion focuses on topography- 
driven transitions on a model surface shaped like a bell 
curve or "Gaussian bump" (see Fig. P), but the same 
mathematical approach can be readily carried over to 
study arbitrary surfaces of revolution that are topolog- 
ically equivalent to the plane. Furthermore, we do not 
expect the results of this analysis to depend qualitatively 
on the azimuthal symmetry of the surface, which is as- 
sumed purely for reasons of mathematical convenience. 

Points on our model surface embedded in three dimen- 
sional Euclidean space are specified by a three dimen- 
sional vector R(r, (p) given by 



where the important a-dependent function l(r) (see Fig. 
|2Jl is defined by 



l{r) = 1 + 



-^5— exp , „ 
Tn \ rf 



(9) 



and it is equal to the radial component of the diagonal 
metric tensor, gap, 



_ 'l{r) 
^"/3 - I r^ 



(10) 



Note that the 5^0 entry is equal to the flat space result 
r^ in polar coordinates while Qrr — lir) is modified in 
a way that depends on a but tends to the plane result 
grr = 1 for small and large r, as illustrated in Fig. |21 

The Gaussian curvature for the bump is readily found 
from the eigenvalues of the second fundamental form [31| , 
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Note that a controls the order of magnitude of G{r) and 
that G{r) changes sign at r = tq (see Fig. ^). The inte- 
grated Gaussian curvature AG'(r) inside a cup of radius 
r centered on the bump is 



AG(r) = 27r 1 - 



1 



(12) 



which vanishes as r ^ 00. Ea. H12|) also shows that the 
positive Gaussian curvature enclosed within the radius rg 
(see Fig. QJ approaches 27r for a ^ 1, half the integrated 
Gaussian curvature of a sphere. 



R(r,, 



r cos( 
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where r and cf) are plane polar coordinates in the xy plane 
of Fig. Q] It is useful to characterize the deviation of the 
bump from a plane in terms of a dimensionless aspect 
ratio 



h 



(7) 



H and t^ 



The two orthogonal tangent vectors t^ 
^ can be normalized to define the Vierbein (orthonor- 
mal basis vectors) E^ and E^ respectively. The com- 
ponents of the spin connection introduced in Eq. Q are 
given by Aq — E^ • 9a E^ J2^ |23|. This leads to a van- 
ishing radial component Ar and 



A, 



(8) 



B. Defect free texture 

For small values of the aspect ratio a, the minimal 
energy texture for the hexatic will be free of defects. The 
ground state configuration 9o{u) satisfies the differential 
equation 



DaD°'9o ~ D'^Aa = , 



(13) 



which results from minimizing the free energy in Eq.Q 
with respect to the field 6{\i) for fixed Aa- When ex- 
pressed in terms of the coordinates in Eq.lO, the solution 
of Eq. (1121 reads: 



9o(u) 



(14) 



where c is an arbitrary constant. The smooth ground 
state texture is thus obtained if the director n forms an 
angle 9o{yi) — —(j)+c with respect to the spatially varying 
basis vector E^ . Note that a solution of the form 9o (u) = 
c represents a defect of charge g = 27r in this " rotating" 
system of coordinates. 
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FIG. 1: (a) The vector field n is confined to a surface shaped as a Gaussian, (b) Top view of (a) showing a schematic 
representation of the positive and negative Gaussian curvature as a background "charge" distribution that switches sign at 
r = ro- Note that, according to the electrostatic analogy, a positive (negative) distribution of Gaussian curvature corresponds 
to negative (positive) topological charge density. 
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FIG. 2: Plot of l( — ) as a function of the dimensionless radial 
coordinate — for a = 1,2,3,4. The arrow is oriented in the 



direction of increasing a. 



As an illustration, consider the projection on the plane 
of the minimal energy texture of an XY model {p ~ 1) 
as shown in Fig. 13 The arrows represent the orientation 
of tilted molecules on this surface in the one Frank con- 
stant approximation. The field clearly displays strong 
frustration along a direction determined by the choice 
of the constant c in Ea. H14|l . If the bump is positioned 
within two very distant walls parallel to the y-axis which 
impose tangential boundary conditions on the molecular 
tilts, the "preferred" direction will be along y jij. The 
texture displayed in Fig. |31can be interpreted as result- 
ing from embryonic pairs of defect dipoles along the line 



X — 0. The distortion energy _Fo of this ground state is 
given by: 



Fo 



1 



Ka / dAg'^'^ido^Oo - Ao.){dp0o - Ap) . (15) 



This expression can be evaluated for an infinitely large 
system by using Eq. I|14|l and the explicit form of the spin 
connection derived in Section [IIAI with the result: 



Fo = ttKa / dr^ j=^ 

Jo i'y'l[r) 



(16) 



It follows from Ea. H16|l that the ground state energy is 
a monotonically increasing function of the aspect ratio, 
proportional to a* for small a. As we shall see, for large 
enough a, it can be energetically preferable to reduce this 
energy by introducing defect pairs into the texture. It is 
convenient to rewrite Ea. (|15|l in terms of the Gaussian 
curvature G{r) and the Green function r(u, u') discussed 
in Appendix 1X1 p^ 



^0- 



Ka 



I dA jdA' G(u) r(u,u') G(u') . 



(17) 



This result is what one obtains by setting all g^ = in 
Eq.Q). The details of the mathematical derivation are 
relegated to Appendix lUl 

Although this result correctly represents the zero tem- 
perature limit of the vector model, corrections may be 
appropriate to describe the physics of ordered phases 
at finite temperature. "Spin-wave" excitations (i.e.. 
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FIG. 3: Projected ground state texture for an XY model on the bump, with the boundary condition that the vector field is 
parallel to the y-axis at infinity. The two insets show the defect pairs suggested by two regions of large frustration, which lie 
close to a circle of radius ro. 



quadratic fluctuations of the order parameter about the 
ground state texture) can be accounted for by integrating 
out the longitudinal fluctuations 0'(u) around the ground 
state conflguration 60 (u) . By letting 9 = 9o + 6' in Eq. Q 



and using Eq.Q we obtain |18ll2l| : 



F = Fo + -Ka 



dAcf^dc^e'dpO' , (18) 



The longitudinal variable 9'{\i) appears only quadrati- 
cally in F and the trace over 9'{vl) can be explicitly per- 
formed with the result 1 321: 



Ve'e^'^ / dA g'-f d^B' d(,0' ^ ^-flFr. 



(19) 



where F^ is the Liouville action, 
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f dA-— f dA fdA' G{u) r(u,u') G(u') . 

(20) 
The first term in this expression is a constant propor- 
tional to the fixed surface area of the frozen topography 
and will be suppressed in what follows. The remaining 
term causes a shift in the coupling constant appearing in 
Eq.dni) from K a to K'^ = Ka- ^^ 'JM^- This "en- 
tropic" correction to the coupling constant Ka at finite 
temperature also arises when defects are present. 

The energy in Ea. p7|) represents an intrinsic, irre- 
ducible energy cost of geometric frustration for textures 



without defects. As we shall see, defects can reduce this 
frustration. However, for small values of a the energy 
cost of this frustration will still be lower than the core 
energies associated with the creation of the unbound de- 
fects and the work necessary to tear them apart. 



C. Energetics of defect pairs on a Gaussian bump 

A quantitative understanding of the energetics of de- 
fects on a fixed topography is essential to calculate the 
critical value(s) of the aspect ratio above which defect- 
unbinding becomes energetically favorable. The first step 
is to calculate the Green's function r(u, u') that governs 
the "coulombic" interaction among defects and between 
each defect and the Gaussian curvature. The inversion 
of the curved space Laplacian can be more easily accom- 
plished by employing a set of "isothermal coordinates", 
such that the resulting Green's function reduces to the 
familiar logarithm of two dimensional electrostatics. As 
shown in Appendix ^ the final result in ternrs of the 
original polar coordinates reads: 



r(u,u') = -^ln[5R(r)2 
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- 23fi(r)5R(r')cos(0-0')]+c , (21) 

where the function 5R(r) can be thought of as a radial co- 
ordinate in the conformal plane resulting from adopting 



an isothermal set of coordinates (see Appendix 
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(22) 



and l{r) is the a-dependent function introduced in 
Eq.(|5J. The constant c depends on the physics at short 
distances, which is discussed in Appendix fXl 

The Green function in Ea. pif) corresponds to free 
boundary conditions at infinity and preserves the cyhn- 
drical symmetry of the metric. It differs from the famiUar 
result in flat space by a non-linear radial stretch corre- 
sponding to a smooth deformation of the bump into a 
flat disk. This Green's function determines an attrac- 
tive interaction for the defect dipole pair. However, the 
Gaussian curvature of the bump also generates a geo- 
metric potential that tries to pull the disclination dipole 
apart. This geometric interaction arises by combining 
cross terms between s(u) and G(u) in Eq.Q with the 
position dependent self-interactions derived in Appendix 
IgI The resulting interaction Fq between defects and the 
Gaussian curvature takes the simple form 

FG=KAj2q.{l~^)V{u,), (23) 

where the geometrical potential V(u) is defined as 



V{u) = - fdAG{u') r(u,u') 



(24) 



The minus sign in front of this geometric potential insures 
that defects of topological charge between zero and An 
are attracted by regions of positive Gaussian curvature 
[30|. For defects with a large topological charge of g > 
47r, the sign of the geometric interaction, Fq, is reversed 
and, and defects of either sign are pushed away from the 
bump. This scenario does not affect the geometry-driven 
defect formation discussed in this paper that relies on 
disclinations whose charge is well below Att. However, as 
a result of the position-dependent self interactions, Fq is 
no longer symmetric under the change q -^ —q, as one 
would expect on the basis of the electrostatic analogy. 
The effect of this asymmetry is small for hexatic order but 
it increase for liquid crystals with p-fold order parameter, 
as p decreases (see Ref. [S^ , and references therein) . 

Gauss' law generalized to curved surfaces (see Ap- 
pendix insures that the geometric force experienced 
by a defect of charge q with radial coordinate r will be de- 
termined only by the net curvature enclosed in a circle of 
radius r centered on the top of the bump. The resulting 
"electric field" is radial as expected from electrostatics 
and is proportional to the gradient of the geometric po- 
tential, which for a Gaussian bump takes the form (see 
Appendix ^ 



yir) = -['^{vw)-i) 



(25) 



For small values of the aspect ratio a the potential in 
Ea. (|25|) can be approximated by 



V{r) 
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(26) 



The resulting force is linear for small r (i.e. near the top 

of the bump) and decays like e ''o for r ^ tq. As the 
aspect ratio a increases, the force generated by the curva- 
ture can overcome the attractive force binding the defect 
pair which varies logarithmically for short distances. As 
a result, oppositely charged defects that were originally 
tightly bound can be separated. 

This argument, however, neglects another complica- 
tion resulting from the curvature of the surface: as the 
aspect ratio is increased, the Green's function r(u, u') in 
Ea. (|21|) and hence the force binding the defects together 
also increases. We illustrate this point in Fig. ^for the 
special case of a positive defect pinned right on top of 
the bump and a negative one free to move downhill at 
r. Upon invoking Equations 10} and 1)21(1 . the potential 
binding the pair, Vpair (r), can be written down exactly 
as: 



^pair V / 



27r \a 



2q^E, 



27r 



dr' 
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(27) 



The first term is the flat space Green's function and we 
have added two disclination core energies. The last term 
represents a "curvature correction" that has the same 
functional form of the geometric potential in Eq. 1)25(1 , but 
it represents a distinct contribution to the total free en- 
ergy. As discussed in Appendix ^ the pair potential 
energy, Vpair [i"), can be understood by applying a gen- 
eralized Gauss' Law to the bump to determine a force 

which is equal to — ^^. The potential follows by inte- 
grating this force along the bump with the length element 
dr^JUr). Although the force is independent of the aspect 
ratio, the length element grows with a (see Fig. 0J, which 
makes the pair more energetically bound for larger values 
of a. A careful calculation of these effects (including the 
contribution associated with the position-dependent self- 
energies) reveals that the geometric force still overcomes 
the binding interaction for sufficiently large values of a. 
An estimate of the critical value of a for which the 
dipole unbinds can be obtained by comparing the mini- 
mal free energy of the smooth frustrated field arising from 
Eg. ((17(1 with the free energy in the presence of defects. 
The latter, as follows from Eq.Q, is composed of three 
contributions: the interactions among the defects, the in- 
teraction between the defects and curvature as given by 
Ea. ((23() and the Gaussian curvature self- interaction. The 
latter is equal to the minimal free energy of the smooth 
frustrated field and is renormalized at finite temperature 
in the same way. (Associated with the cutoff is a micro- 
scopic core energy, E^, that we expect to be independent 
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FIG. 4: EfTect of changing the aspect ratio of the bump on 
the work needed to puU apart two oppositely charged defects. 
Positive and negative defects are represented by open circles 
with their sign printed. The line elements corresponding to 
the projected length dr for the two aspect ratios are shown. 



FIG. 5: Geometric potential V{r/ro) as a function of the 
dimensionless ratio of r and ro for q = 1, 2, 3, 4, 5. Note that | 
V{r) |<^ I V(0) I for r > rg. The arrow points in the direction 
of increasing a. 



of the defect position on the bump, as long as the radius 
of curvature is much greater than any microscopic length 
scale.) The remaining two contributions are not renor- 
malized by thermally induced spin wave fluctuations [2lj . 
As shown in Appendix El the difference in free energy of 
a defected texture described by Eq.l^J relative to the 
defect- free result Eq. ljlTI) can be written as: 

Ar (a) 1 ■^-^ ■^-^ 

~^ ^ ol^l^ ^^J-rair,, 0j, Tj, 0j) 



Ka 



1=1 j^i 



Na j^ Na 



(28) 



i=l 



where we have assumed overall charge neutrality for the 
defect configuration. The subscript in F^ indicates that a 
constant microscopic core radius a has been absorbed in 
the definition of the Green function so that the argument 
of the logarithm in Ea. l|21|) becomes dimensionless, as in 
Eq.dnSJ, 



Ta{ri,(t>i,rj,(l)j) =r(r,;,(/),,rj,0j) + — Ina . (29) 

This microscopic cutoff a corresponds to a constant core 
radius for each defect and it is of the order of the spacing 
between the microscopic degrees of freedom. The sum of 
microscopic core energies in the fourth term of Ea. H28|) 
needs to be fixed phenomenologically or from models that 
go beyond simple elasticity theory j^. Note that both 
'^a{'i'i,4>i,rj,4>j) and V{r) depend on a. If a becomes 
sufficiently large so that AF is less than zero, one or 
more disclination dipoles unbind in the hexatic phase at 
a sequence of critical values a^ ■ The analogous defects 
in XY-model textures of tilted liquid crystal molecules 
would be +/— vortex pairs. 



III. CURVATURE INDUCED DEFECT 
GENERATION 

A. Onset of the defect-dipole instability 

If a dipole is created, say, along the line r = ro of 
zero Gaussian curvature, the positive disclination will 
be pulled towards the center by the positive curvature 
while the negative one will be repelled into the region of 
negative curvature. The net result is a reduction of the 
total free energy of the order of the depth of the poten- 
tial well since the logarithmic binding energy is approx- 
imately constant compared to the geometric potential. 
An approximate analytical treatment is obtained by as- 
suming that the positive defect sits right at the center of 
the bump and the negative one at a distance of the order 
of ro- The validity of this approximation scheme can be 
checked by numerically minimizing the energy with re- 
spect to the position of the defects, as discussed in Sec- 
We assume charge neutrality so that the 



tion ljIirB|l 

two defects have equal and opposite topological charges 
of magnitude q. For order parameters with a p-fold sym- 
metry, the minimal topological charge is q = -^. The 
approximate free energy cost to generate this defect pair 
then follows from Equations l(^ . lt?7|l and (J2H1): 



AF(a) 
Ka 



2tt 



1 



ln(^ 
a 

q 



V{r) 



i^J^W 



q[l 
^ Ka 



i^)"'") 



(30) 



The internal consistency of the formalism can be checked 
by investigating the limit r —^ a. As the negative defect 
approaches the positive one at the center of the bump, we 
have V{a) « V{0) and the energy tends to the expected 
flat space result 2q^Ec. 

The equilibrium position of the negative defect turns 
out to be for r equal to a few tq, so the terms containing 
V{r) in Ea. H3U|) can be dropped as a first approximation 
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FIG. 6: Plot of the geometric potential evaluated at the center 
of the bump, V{0), for aspect ratios a between and 5. The 
continuous line is plotted using the exact form of V{r) while 
the dashed line is obtained from the low a expansion given in 
Eq.lUni- 



because V{r) decays exponentially (see Fig. IS}: 



AF{a) 
Ka 



2tt \a 



9 1 



4^J^(«) 



E, 



2,'^^AU) 



To estimate the critical value of the aspect ratio for which 
the first dipole unbinds, let r '^ ro and solve for cuc in 
Ea. (|31|l . Because a different choice for the core energy 
Ec can be accounted for by rescaling the core size a in 
Ea. (|31|l . the condition for unbinding is 



\vm > 



2q 



{An - q) \a' 



with 



(32) 



(33) 



If we know -jr^ and ^, the critical aspect ratio ac can 
be obtained from inspection of Fig. Elwhere V^(0) is plot- 
ted as a function of a. The Taylor expansion of V{r) 
derived in Eq.(|2SJ| gives V{0) « — ^ to leading order 
in a. Inspection of Fig. shows that this approxima- 
tion works sufficiently well even for aspect ratios of order 
unity. Upon substituting for V{Q) into Ea. H32|l . we obtain 
an estimate of how Uc depends on ^ : 



{An - q) \a' 



ln(^ 



(47r - q) 



AttE, 



Ka 



(34) 



Note that a critical height he = aero for defect unbinding 
is predicted for fixed ^ with defect charge g = — for 
all integer values of p. The validity of this approximate 
relation is tested in Section IlIIBI 

The continuum theory adopted here is valid in the limit 
ro ^ a. If Ec can be neglected compared to Ka, the de- 
fect unbinding instability is triggered when the energy 



gain derived from letting the defects screen the Gaus- 
sian curvature (approximately given by ^^^(O)) overcomes 
the work needed to pull them apart a distance tq. This 
work, of order fj^ In (^) , increases very slowly with large 
— , hence the continuum approximation can be satisfied 
while keeping the work finite. Note that the result does 
not depend on the size of the system R because we as- 
sume overall disclination charge neutrality and the as- 
sumption that R ^ tq. In this limit, boundary effects 
can be ignored provided that they do not impose a topo- 
logical constraint on the phase of the order parameter. 
An aligning outer wall in a circular hexatic sample, for 
example, would force the bond angle field to rotate by 
27r, leading to six defects in the ground state even in flat 
space. The interesting physics which results is addressed 
in Section IW1 



B. Numerical investigation of defect-unbinding 
transitions 

The disclination unbinding transitions can be inves- 
tigated more quantitatively by minimizing numerically 
AF{a) in Ea. pSfl with respect to the positions of the de- 
fects. The aspect ratio above which AF{a) becomes neg- 
ative corresponds to the threshold value ac (analogous to 
a first order transition) for which the singular field is en- 
ergetically favored with respect to the smooth texture of 
Fig. 131 We emphasize that the energy landscape can 
have two minima. The first occurs when two oppositely 
charged defects form a closely bound dipole (with sepa- 
ration of the order of the cutoff a) and hence annihilate 
each other leaving a smooth texture. The second min- 
imum corresponds to an unbound pair (with separation 
of a few ro) and it disappears when the geometric force is 
too weak to overcome the binding force of the pair. This 
scenario occurs for a finite value of a characteristic of 
the geometry of the substrate above which the formation 
of an unbound dipole is possible, albeit energetically un- 
favored (see Fig. Cj). As a is increased above ac, the 
smooth-texture minimum becomes metastable and the 
unbinding of a defect pair is the most likely scenario. 

It is useful to parameterize AF(a) in terms of the di- 
mensionless radial coordinates f,- = — . The geomet- 
ric potential is defined in Ea. (jA9(l as a function of fi, 



V{r) 



Vai^^,)^ where 



Vcix) 



y 



(Vl + aVexp(-y2)-l) . (35) 



In order to write the defect-defect interaction in terms 
of the dimensionless radial coordinate fi , we introduce a 
new function ^{fi) defined by 



sR(n) = — exp[y(ri)] = — exp 
ro ro 



Vo, 



SR(r.) 
ro 



(36) 
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FIG. 7: Plot of /S.F/ Ka versus x\ and X2 the positions of the 
negative and positive defects respectively in units of ro, for 

a = 1.2. A constant energy offset equal to — fir^'^lo/) '^^^ 
been neglected. The metastable minimum at xi ~ 1.3 ro and 
X2 — 0.2 ro corresponds to an unbound pair (see Fig. |HK)- As 
a is decreased further the energy barrier that separates this 
minimum from the smooth texture solution (corresponding to 
xi approaching X2) disappears and the two opposite defects 
annihilate. 



where Ea. ljA8|l was used in the last step. We can now 



transform ra{ri, c 



j) by ehminating 5R(ri) in favor 



of ^(ri). Thus we have using Equations (|C21|I and (jA12p 






Nd 



I : r« : ' ]tVJ ) 



-^(lt'lj^{n,(l)i,rj,(l)j) 
±Y:ifH-), (37) 



2n 



where we have exploited charge neutrahty and Ea. (|C22|l . 
The free energy minimized with respect to the positions 



of the defects, min 
to Eq.lEHI), as 



AF 



Ka 



can now be written, according 



Ka 






where 



/(«) 



mm 



A'd 






Jv^* 



i:*(i-ii)n. 



(39) 



Note that in the second term in Ea. (|38|l the core energy 
of each defect has been absorbed in the modified core ra- 
dius a', as defined in Ea. (|33|l . This generates an energy 




FIG. 8: The equilibrium defects positions are illustrated 
schematically in the case of one (a) and two dipoles (b) . We 
assume free boundary conditions at infinity, as in Fig. |3| so 
that the effect of image charges can be neglected. 




FIG. 9: Plot of /(q) versus the aspect ratio a obtained by 
minimizing over the single-dipole defect configuration repre- 
sented schematically in Fig. |H^. As discussed in the text, 
the first unbinding transition occurs when / (Qci ) is equal to 
— 1^ In (^) . The value Uci is indicated by the dashed line for 
^ = 10* and q = ^ . Note that no minimum (corresponding 
to an unbound pair) exists for a less than 1 (approximately) , 
hence the curve cannot be continued to the origin (see Fig|7J. 



cost for unbinding that can be overcome if f{a) assumes 
sufficiently large negative values. Hence it is sufficient 
to study numerically how AF{a) varies as a function of 
a single parameter, eg. ^. As an illustration of this 
approach, we study explicitly the unbinding of one and 
two disclinations pairs leading to the ground states repre- 
sented schematically in Fig. |H1 The smooth ground state 
becomes unstable to the formation of one defect dipole 
first. The critical aspect ratio above which this scenario 
occurs can be determined with the aid of Fig. |^ where 
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FIG. 10: Plot of /(a) versus the aspect ratio a obtained 
by minimizing the two-dipole defect configuration represented 
schematicaUy in Fig. Et . The second unbinding transition 
occurs when f{ac2) is equal to — ^In (-^)- The value ac2 
is indicated by the dashed line for -^ = 10 and q — ^ 
(compare with Fig. |^. 



the function /(a) introduced in Ea. (|39|l is plotted as a 
function of the aspect ratio. From Ea. H38|l we see that 
Qfci is determined by 



/KJ--|^ln(5 



(40) 



For. 



10"* and q = ^, we obtain a critical aspect ratio 
ttci ~ 3.2. As a comparison, the approximate condition 
derived in Eg. 1)32(1 gives ac^ « 3 when used in conjunc- 
tion with Fig. El The rougher estimate in Ea. (|34|l leads 
(for q = ^) to aci ~ 2.6. This discrepancy is easily un- 
derstood considering that Eq. ((34|l was derived by means 
of a low a expansion. 

The critical aspect ratio ctci is too low for the two- 
dipole defect configuration to become energetically fa- 
vorable with respect to the smooth ground state. Indeed, 
inspection of Fig. ^| reveals that the critical aspect ra- 
tio ac2 for which the "two-dipole instability" sets in is 
approximately equal to 3.6 for the same choice of pa- 
rameters used in the single pair case. Note that, in the 
presence of two dipoles, the energy cost arising from the 
second term in Eq. ll5S|l is twice as large because there are 
four defects rather than two. However, for a > 4.2 gen- 
erating two dipoles becomes more energetically favored 
than a single dipole (see Fig. [TT|) . The approach illus- 
trated here can be used to calculate a cascade of defect 
unbinding instabilities at critical aspect ratios a^ involv- 
ing higher number of dipoles and their equilibrium con- 
figurations in the ground state. Note that the unbinding 
eventually stops since the integrated Gaussian curvature 
in the top cup cannot exceed 27r. We expect the quali- 
tative features of our analysis to be independent of the 
exact shape of the bumpy substrate although the specific 
values Uci depend on the geometry and the choice of the 
microscopic parameters a and E^- Finally, we empha- 
size that the curvature-induced unbinding is similar to a 



FIG. 11: Plot of ^- "^ versus a corresponding to a single 
dipole (continuous line) and two dipoles (dotted line) for ^ = 
10*. The critical aspect ratios a^ and Oca are indicated by 
dashed lines. Note that the aspect ratio for which the two 
dipole configuration becomes energetically favored occurs for 
a > 4.2. 

first order transition and occurs for rather pronounced 
deviations from flatness (ie. large a). 

C. Single vortex instability 

The unbinding of defect pairs may not be the most 
likely scenario if the size of the system R is sufficiently 
small. In this case, the creation of a single vortex at the 
center of the bump may become energetically favorable 
for lower aspect ratios than required by the defect dipole 
instability. The equation for the bond angle field 9s(u) 
for a single defect of charge q at the center of the bump 
is given by: 



m 



_q_ 

2tt 



-1 



(41) 



where the bond angle is measured with respect to the 
rotating basis vectors corresponding to the polar coordi- 
nates discussed in Section III Bl Upon substituting 6's((/)) 
in Eq. Q and subtracting the free energy Tq correspond- 
ing to the defect-free texture we obtain: 



AF(a) 
Ka 



in 



In 



5R(i?) 



9 1 



47r 



1/(0) 



, Ec 



(42) 



where Ec was added by hand. The same result is ob- 
tained by using the more general formalism developed in 
Appendix D. Indeed, by letting the position of an iso- 
lated defect tend to the center of the bump in Eq. (|D24|) 
we obtain the energy of the singular field in the case of 
free boundary conditions and the result matches Ea. (|42|l . 
As discussed in Appendix ^ a defect located at r^ is 
attracted to the boundary at R for free boundary con- 
ditions. One can think of this interaction as resulting 
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from an image defect of opposite sign behind the edge of 
the sample at position r'^ such that the following relation 
holds in terms of the conformal radius 3? (r') 



W^) 



^{Rf 



(43) 



This result can be understood by analogy to the familiar 
electrostatic problem of a charged line located a distance 
Ti from the center of a cylindrical grounded conductor 
whose axis is parallel to it [SJI- The analogy becomes 
precise if one lets r.i — *■ '^{vi) as explained in Appendix 

El 

If the geometric potential is not strong enough (as in 
the flat space limit a = 0), the defect will migrate to the 
edge of the sample and annihilate with its image leaving 
a smooth field. On the other hand, when the aspect ra- 
tio is sufficiently large, the defect can lower its energy by 
sitting at the center of the bump. Comparison of Ea. l|I^ 
with Ea. H31|l shows that, unless R ^ r^, the energy of 
the single vortex instability will be lower or at least com- 
parable to the unbinding of a defect dipole. In fact, the 
threshold as that a needs to exceed to trigger the single 
defect instability is easily obtained if the values of the ge- 
ometric potential at the origin are tabulated for different 
aspect ratios, as illustrated in Fig. [S] The condition for 
single vortex generation reads 



\vm > 



{q - 47r) 



ln(- 



(44) 



Using the same method adopted to derive Ea. H34|l we 
obtain an estimate of how as depends on -^ (compare 
with Eq.lISlI)): 



{4n-q) \a' 



(45) 



The single vortex instability is reminiscent of vortex gen- 
eration in rotating superfluid helium with a playing the 
role of the angular speed fl. For a volume of helium con- 
tained in a cylindrical vessel of radius R and rotating 
uniformly with constant angular speed, the critical value 
r^ci above which defect generation occurs is given by [33 



nr 



K (R 



(46) 



where K — ^2iL ig the magnitude of the quantum of 



circulation and a the core radius 143 



Note that Q.c 



decreases as R increases, unlike as which diverges log- 
arithmically. Thus, the single defect instability studied 
here is a finite size effect. In contrast, the disclination 
unbinding studied earlier in this section does not depend 
on the system size because of charge neutrality. Hence 
the thermodynamic limit can be safely taken, provided 
the characteristic length over which the curvature varies 
(ie. ro) is not too large compared to a (see Ea. (|5Hl '). 

In considering the case of small system size, it is im- 
portant to keep in mind two assumptions implicit in the 



present treatment. The radius of curvature — must be 
much larger than the core radius everywhere for the con- 
tinuum approach to be valid, that is tq 3> a a. Addition- 
ally, the Gaussian curvature must be vanishing small at 
the edge of the system which requires R to be larger than 
a few tq. 



D. Lattice of bumps, valleys and saddle points 

In some experimental realizations perhaps modelled 
on those of Ref. [ij the topography will be periodic. 
In this Section we discuss qualitatively how the results 
described above generalize to a 2-dimensional lattice of 
bumps with variable aspect ratio for both square and 
triangular lattices. A more quantitative approach to this 
problem would involve finding conformal set of coordi- 
nates for periodic boundary conditions. This is possible 
in principle but more involved since cylindrical symmetry 
is now lost. Nonetheless, the intuition gained by study- 
ing the single bump allows us to make some guesses for 
the ground state. We first note that the geometric poten- 
tial generated by the lattice of bumps is not simply the 
superposition of results for single bump potentials. This 
is caused by the non linear relation between the surface 
height and the Gaussian curvature acting as a source for 
the geometric potential. To explore this point further, 
consider what happens when four bumps are placed at 
the vertices of a square. At the center of the square a 
minimum of the height function occurs corresponding to 
a new region of positive Gaussian curvature. This effect 
is particularly acute for rg < L, where L is the bump 
spacing. In general interference between bumps creates 
a dual lattice of valleys. A similar breakdown of the su- 
perposition principle arises for triangular lattices. 

As the aspect ratio of hilly landscapes such as those 
shown in Fig. E| is increased, defects can be created 
to screen the Gaussian curvature. Their positions can 
be guessed by considering a unit cell of the lattice such 
that the integrated Gaussian curvature vanishes. For a 
square lattice, we conjecture that the first topography 
induced transition is associated with the appearance of 
positive defects at the top of the bumps and negative 
ones half way between them in the vertical or horizon- 
tal direction (see Fig. [T^ . This two- fold degeneracy is 
compatible with the symmetry of the lattice and anal- 
ogous to the freedom in choosing the axis along which 
the first disclination-dipole appears on the single bump. 
The negative defects are shared between two adjacent 
cells while the positive ones are shared among four cells 
thus ensuring overall charge neutrality. As the value of 
a increases even more, one might expect an additional 
positive defect appears in the valley located at the center 
of each cell and two additional negative defects shared 
with the adjacent cells are created between the bumps 
at right angles to the direction discussed above (see Fig. 

For the triangular lattice, we conjecture that the first 
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FIG. 12: (Top) Ground states for square (left) and triangular 
(right) arrays of bumps. The first and second rows correspond 
to moderate values of the aspect ratio a respectively. For 
simplicity, we assume that ro, the bump width is comparable 
to the lattice spacing. Positive defects (red dots) "screen" 
regions of positive Gaussian curvature while negative ones 
(blue dots) are located on the saddles of the "hilly" landscape. 



transition corresponds to positive defects on top of the 
bumps and negative ones between the bumps along one of 
the three axis of symmetry of the unit cell. As the value 
of the aspect ratio is increased, additional positive defects 
appear on the six minima of the surface and negative ones 
are generated along the remaining two axes of symmetry 
of the unit cell (see Fig. I12|l . A simple count of the total 
defect charges enclosed in the unit cell shows that this 
scenario also satisfies the requirement of defect charge 
neutrality. 



IV. 



DEFECT DECONFINEMENT 



With potential experiments in mind [23J , it is interest- 
ing to consider the case of hexatic order on a bump encir- 
cled by a circular wall of radius i? 3> rg which aligns the 
hexatic bond angles. As a simple model, imagine an array 
of hexagons which locally achieve a common orientation 
tangential to the wall (see Fig. 118b). The hexatic order 
parameter will thus rotate by 27r upon making a circuit of 
the wall, insuring that at least six defects of "charge" ^ 
must be included in the ground state for all values of the 



aspect ratio. These boundary-condition induced defects 
will interact with the Gaussian curvature of the bump 
and with the wall. The Nd defects contribute large (con- 
stant) self energies of the form X]j=i ^Aqf In (■§ ) that 
dominate the total energy for sufficiently large systems. 
Since X)i=i 9* must be equal to 27r, the energy is mini- 
mized when the defects split up into the smallest possible 
charges. 

The equilibrium defect configuration must minimize 
the free energy taking into account the confining po- 
tential generated by the Gaussian curvature and the in- 
teractions of the defects with the boundary and among 
themselves. The repulsive force exercised by the wall 
on a defect located at 5R(ri) in the conformal plane can 
be computed by placing an image defect of same charge 
outside the wall at position ,S, '^ . The mathematics re- 
sembles the problem of finding the magnetic field of a 
line current located at a given distance r.i from the cen- 
ter of a cylinder of high permeability material and whose 
radius R is greater than r^ ^J . The analogy is complete 
upon performing the change of coordinates r — > 3?(r) 
and identifying the gradient of the bond angle i9q0(u) 
with the magnetic field. This is explained in detail in 
Appendix O and IdI where we introduce a conjugate func- 
tion x(u) analogous to the vector potential that simplifies 
the analysis of this problem. Thus, each of the Nd defects 
will also interact with an equal number of image defects. 
This situation can be described mathematically by de- 
riving an appropriate Green's function F^ that includes 
the images, as discussed in Appendix IdI (see Ea. HD19|) ). 
The resulting free energy F^ reads: 



pN ^ ^d 
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(47) 



where Fq is defined in Ea. l|17|l and the Green's function 
r^{xi]Xj) is given by 



i yXi^ Xj ) 



1 

47r 



In [xj + X, — 2xiXj cos {(j)^ — (j)j)] 



in 



In [x'^x'j + 1 — 2xiXj cos {(f>i — (j)j)\ 



(48) 



The last term accounts for the interaction with the im- 
age defects and the superscript N indicates Neumann 
boundary conditions on an appropriate potential func- 
tion. Here, we use scaled coordinates in the conformal 
plane Xi = wifl • The interaction of the defects with the 
curvature is not affected by the presence of the distant 
wall. 
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FIG. 13: (a) The free energy for a nematic (double headed 
vector field) living on a Gaussian bump surrounded by an 
aligning circular wall is plotted for a = 2 as a function of the 
scaled radial coordinates xi and X2 of the two disclinations. 
The radial coordinates have been scaled by ro and the size of 
the system is _R = 7ro. Note that the energy plot is symmetric 
with respect to the line xi = X2- (b) Schematic illustration 
of the positions of the two disclinations (black dots) corre- 
sponding to the deep energy minima at positions xi — 0.04 
and X2 = 4.9 (or viceversa) and to a shallow minimum at 
xi = 2:2 = 4.7. The two defects are on opposite sides of 
the bump. The continuous line corresponds to the circular 
boundary while the dashed one to the circle of zero Gaussian 
curvature and radius ro (drawing not to scale). 



To provide an illustration of the combined effect of 
curvature and boundary conditions on tangential vector 
order, we first consider the simpler case of a nematic 
order parameter with periodicity equal to tt. This sim- 
plified model neglects differences in the elastic constants 
for bend and splay and does not incorporate any effect 
due to the uniaxial coupling of the nematogens to the 
curvature. In this case, minimization of the logarithmi- 
cally diverging part of the free energy (fourth term in Eq. 
I47|l suggests that there will be only two disclinations of 
charge q = it displaced along a radial direction (see Fig. 
113b). By applying Ea. H47|l . we can parameterize the en- 
ergy of the system in terms of the scaled radial coordi- 
nates, xi and X2 of the two disclinations. The resulting 
energy landscape is plotted in Fig. E| (for a = 2 and 
R — 7 To) and clearly reveals two minimal-energy config- 
urations. The first minimum corresponds to one discli- 
nation confined at the top of the bump (slightly shifted 
from the center) and the other at a radial distance ap- 
proximately 70% of R (see Fig. 113b bottom panel) . The 
second minimum corresponds to a fully deconfined state 
with both disclinations placed symmetrically at approxi- 
mately 67% of R (see Fig. [T^ top panel). As the aspect 
ratio is raised even further, the saddle in the energy land- 
scape of Fig. 113b becomes a minimum corresponding to a 
configuration in which both disclinations are confined in 
the cup of positive Gaussian curvature by the geometric 
potential. 

As illustrated in Fig. ^1 there is a critical value of the 
aspect ratio, ao — 1.5, above which it is energetically fa- 
vorable for the system to have one disclination confined 
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FIG. 14: Plot of the free energy of a nematic (double headed 
vector field) on a Gaussian bump encircled by an aligning wall 
as a function of a. The dotted line represents the energy of the 
fully deconfined configuration in Fig. 113b top panel while the 
continuous line corresponds to the defect pattern illustrated 
in the bottom panel of Fig. 113b . The energy of the fully 
deconfined configuration is approximately independent of a 
because the two disclinations are far away from the bump. 



at the top of the bump. For a < ao the fully decon- 
fined configuration becomes energetically favorable, but 
the two minima can still coexist. As a is decreased even 
further, the repulsion between the two disclinations over- 
comes the confining force of the geometric potential and 
makes the second minimum in Fig. 114b (corresponding 
to the partially confined configuration) disappear alto- 
gether. This "spinodal point" occurs for a « 0.9 on the 
Gaussian bump. The specific values of the critical aspect 
ratios are geometry dependent, but the generic mecha- 
nism of deconfinement depends only on a large separation 
of the length scales rg and R that control the interaction 
with the curvature and the boundary respectively. 

The analysis for the hexatic case is complicated by the 
fact that more defect configurations are possible when 
six defects are present. We start by noting that even 
in flat space (a = 0) there are two natural low energy 
defect configurations with high symmetry: the ground 
state corresponding to the six defects sitting at the ver- 
tices of an hexagon and a higher energy state given by 
a pentagonal distribution of defects with the sixth de- 
fect sitting at the center of the circular sample (see Fig. 
115b). As the aspect ratio is raised, the pentagonal ar- 
rangement becomes energetically favored since it pays to 
have a defect confined in the (geometric) potential well 
at the origin (see Fig. llSb). To study the transition, it 
is useful to derive expressions for the energy of the two 
defect configurations as a function of the radius of the 
outer defect ring r. Every defect (except the one at the 
origin, possibly) has the same scaled coordinate Xi — x 
and the angles between two defects are integer multiples 
of — where n is the number of defects in the outer ring 
{n — 5 for the pentagon and n = 6 for the hexagon). In 
this case, the sums involved in the first (interaction) term 
in Eq. 147|) can be efficiently evaluated using the following 
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identity: 



_. n— 1 

2^ 



In 



27ri 

1 — 2pcos( ) 

n 



= ln(l-p")-ln(l-p) 



(49) 



Upon using Ea. (|49|l with p = I and p ^ x^ to evaluate 
the sums arising from the first and the second term of 
the Green's function in Ea. H48() respectively, we obtain 
the free energy Fh for the hexagonal configuration: 



■Fg(a) _ _7r 
Ka " 6 '' 



/jR(r)y /jR(r)\ 

U(i?); Uw; 
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IItT , , TT 



5R(i?) 



(50) 



The free energy for the pentagonal configuration Fp is 
readily obtained after similar manipulations 
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FIG. 15: Plot of the free energy of an hexatic phase (draped 
on the Gaussian bump encircled by a wall) as a function of 
a. The dotted line represents the energy of the hexagonal 
configuration illustrated in the top panel on the left while the 
continuous line corresponds to the pentagonal arrangement in 
the bottom panel. The outer defect rings in both configura- 
tions are approximately 90% of R. The critical aspect ratio 
od corresponding to the deconfinement transition discussed 
in the text is indicated by the dashed line. 
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^{R) 



.(51) 



Note that these manipulations are very similar to the 
ones necessary to describe superfluid helium in a cylin- 
der of radius R [2a|- In fact, the superfluid problem is 
analogous to the case of hexatic order with free bound- 
ary conditions on a circular boundary of radius R (see 
Appendix ^J . The rather unusual form of the argument 
of the logarithm in Equations 1)50(1 and 1)51(1 arises from 
the sum over the image defects whose positions depend 
non-linearly on the position of the defects themselves. 

Minimization of Equations 1(50(1 and ((51(1 with respect 
to r fixes the distance of the outer defects. The resulting 
minimal energies Fp and Fh are plotted as a function of 
a in Fig. 115b .. The critical value an , for which Fp < Fh 
can be easily estimated by realizing that Fh is approxi- 
mately independent of a because the disclinations are far 
from the bump. On the other end, Fp decreases with a 
because the confined disclination is trapped in a poten- 
tial well whose depth is approximately given by — j^ 7ra^ 
(see the second term of Eq. ((5111 and the low a expansion 
for V{Q) derived in Eg. ((26(1 ). The critical aspect ratio, 
ao, for which the deconfinement transition occurs can 
be estimated by setting the depth of this potential well 
equal to the energy difference between the hexagon and 
pentagon configurations in flat space. The latter can be 
read off from the energy diagram in Fig. 115b and the 
result is approximately 0.3Ka which leads ao ^ 1.1 in 
agreement with the value indicated in Fig. 115b . 

As the aspect ratio is raised even further, less sym- 
metric defect configurations become energetically favored 
corresponding to a larger number of disclinations con- 
fined in the cup of positive Gaussian curvature. For ex- 
ample, when two disclinations are confined within r = tq, 



the outer defect ring is given by four defects approxi- 
mately located at the vertices of a square. We note that 
these defect configurations cease to exist at low aspect 
ratios because they require the geometric potential to 
overcome the strong repulsive interaction between the 
confined defects. As discussed earlier for Fig. lfni) . the 
actual values of the aspect ratios involved depend on the 
specific geometry of the substrate. However the basic 
mechanism behind the deconfinement transition is more 
general. 

Note that, as a increases, the geometric mechanism of 
defect dipole unbinding discussed in the last section may 
also set in. Because of the presence of one or more pos- 
itive defects at the top of the bump, the critical aspect 
ratio necessary to unbind one dipole will be larger than 
what was calculated before. If dipole unbinding does oc- 
cur, the new defects will "decorate" the existing patterns 
by adding new positively charged disclinations in the re- 
gion of positive Gaussian curvature and expelling the 
negative ones in the external region of the bump (r > tq) 
where the Gaussian curvature is also negative (see Fig^l. 



V. CONCLUSION 

We have discussed how the varying curvature of a sur- 
face such as a "Gaussian bump" can trigger the genera- 
tion of single defects or the unbinding of dipoles, even if 
no topological constraints or entropic arguments require 
their presence. This mechanism is independent of tem- 
perature if the system is kept well below its Kosterlitz 
Thouless transition temperature. It would be interest- 
ing to revisit Kosterlitz-Thouless defect unbinding tran- 
sitions on surfaces of varying Gaussian curvature in the 
presence of a quenched topography [I3 in the light of the 
present work. One might also explore the dynamics of 
the delocalization transition that occurs when a bump is 
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confined by a circular edge and the aspect ratio is lowered 
until the defects, initially confined on top of the bump by 
the geometric potential, are forced to "slide" towards the 
boundary. Quantitative studies of periodic arrangements 
of bumps would be interesting and could be inspired by 
fruitful analogies with methods and ideas from solid state 
physics. 

We also hope to extend this work by considering crys- 
talline order on bumpy topographies and taking ex- 
plicitly into account the screening of clouds of disloca- 
tions and possible generation of grain boundaries |16J| . 
Such an analysis would facilitate comparison with exper- 
iments performed with a single grain of block copolymer 
spherical domains |4J| on a suitably patterned substrate 




FIG. 16; Graphic construction of the stereographic projec- 
tion. Regions close to the north pole have larger images in 
the conformal plane than regions of equal areas close to the 
south pole. The stereographic projection preserves the topol- 
ogy of the surface provided all points at infinity are identified 
with the north pole. 
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APPENDIX A: GREEN'S FUNCTION AND 
ISOTHERMAL COORDINATES 

The analysis of ordered phases on curved substrates 
can be simplified by rewriting the original metric of the 
surface gabi'^) in terms of a convenient set of coordinates 
x(u) = (a;(u),y(u)) such that the new metric gab{x) 
reads 



i(x) = eP'^^^y^Sa 



(Al) 



The metric gab differs from the flat space one Sab only by 
a conformal factor e^'^^'^^ that embodies information on 
the curvature of the surface 31] . These isothermal coor- 
dinates can be used to map arbitrary corrugated surfaces 
onto the plane ^7!|. The mapping is conformal so angles 
are left unchanged but areas are stretched according to 
the position-dependent conformal factor e''^^'^^. A famil- 
iar example is provided by the stereographic projection 
that maps a sphere onto the conformal plane as illus- 
trated in Fig. ^1 The Green's function assumes a very 
simple form after a conformal transformation, because 
the Laplace operator reduces to the familiar flat space 
result when expressed in terms of isothermal coordinates. 
In what follows, we demonstrate that this transformation 
provides the basis for an efflcient strategy to determine 
the Green's function on a bumpy substrate. We start by 
deriving the radial change of coordinates 3fi(r) that trans- 
forms the original metric of the Gaussian bump, i.e.. 



ds^ = 1 



dr^ 



(A2) 



into the locally flat metric (in polar coordinates), 

where p(r) and 3?(r) are independent of the azimuthal 
coordinate (j) because of cylindrical symmetry. This met- 
ric is equivalent to gabi^) upon switching from cartesian 
(x,y) to polar coordinates (3?(r), (j>). To simplify the nota- 
tion we introduce the a-dependent function l{r) defined 
by 



l{r) 



1 



■ exp 



(A4) 



and plotted in Fig. [21 for different choices of a. 

The equivalence of the metrics in Equations (jA2|) and 
(IA3(I requires that 3fi(r) satisfies the differential equation 



(A5) 



dn V 


-——dr . 
r 


The conformal factor is thus 


given by 


gP(r) ^ 


(s)^^ 


The solution of Eq.(|A5|| is 




^{r) ^ Are-^r 


= 1s^(V'('-') -1) 



(A6) 



(A7) 



where it is convenient to set the arbitrary constants A 
and c to unity and infinity respectively. This non-linear 
stretch of the radial coordinate leaves the origin and the 
point at infinity invariant and can be concisely written 
as 



sR(r)=re^('') 



where the function V{r) defined by 



Vir) EE - 



dr 



(A8) 



(A9) 
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plays an important role in our formalism and its inter- 
pretation as a sort of geometric potential is explored in 
detail in Appendix IbI 

The Poisson equation for the Green's function r(u, u') 
on a surface with metric tensor gap and point source 
S{u,u') reads ^: 



D"DaT{u,u') 



^/5 ' 



(AlO) 



where the covariant Laplacian is given for general coor- 
dinates by: 

D'^D^ EE (1/V^)9„[V55"^5^] • (All) 



The conformal change of coordinates transforms \/g{r, (f) 
into 6^^^3(1701 and g°'l^ir,(j)) into e-PW5"'3(3fi, </)). 
The factors of e''*^''^ inside the square brackets in Ea. ljAll|l 
then cancel and we are left with the flat space Lapla- 
cian in the polar coordinates (5ft(r ),</)). We conclude 
that r(u, u') is simply the Green's function of an unde- 
formed plane expressed in terms of the polar coordinates 
(5R(r),0): 



r(u,u') 



1 

47r 



ln[^{ry+^{r'y-2n{r)n{r')cos{(j>-4>')] 



(A12) 

where an arbitrary additive constant C (which can be 
used to satisfy boundary conditions at infinity) has been 
dropped. We note that r(u, u') differs from the flat space 
Green's function by a non linear stretch of the radial co- 
ordinate. In Appendix El we will use r(u, u') to solve 
Poisson's equation on an infinite bumpy domain and cal- 
culate the energy stored in the field. As in flat space, the 
Green's function r(u, u') will be suitably modified in a 
finite system in a way that depends on the boundary con- 
ditions chosen at the edge of the sample (see Appendix 

We conclude this appendix by evaluating r(u, u') when 
the two points u and u' are assumed to be separated by a 
fixed distance a small enough so that the surface can be 
approximated by the local tangent plane to the Gaussian 
bump. This short distance behavior will be useful when 
evaluating the effect of a constant core radius on the 
energetics of a disclination at an arbitrary position. The 
fixed microscopic length a on the bump is stretched when 
projected in the conformal plane (see, e.g., Fig. I16|) and 
assumes the position dependent value X{x,y) given by 



A(x,y) 



(A13) 



For a Gaussian bump cylindrical symmetry requires that 
X(r, (j)) is dependent only on r and can be explicitly writ- 
ten upon using Equations ljA6|) and IjASp as 



A(r,(/.) 



^{r) 



= ae^W. 



(A14) 



We can now evaluate r(u, u') in the limit u' -^ u -I- a, 
where this concise notation means that the two points on 



the surface with coordinates u and u' are separated by an 
infinitesimal distance a measured on the bump. It does 
not matter in what direction the two points approach 
each other as long as a is small compared to the local radii 
of curvature. Upon using Equations (|A12|1 and (|A14|) . we 
obtain 



r(u,u + a)^-^ln(a^)-^ 



(A15) 



We note that r(u, u + a) for fixed a is not a constant 
like in fiat space but varies with position as the function 
V{r), reflecting the lack of translational invariance on 
an inhomogeneous surface, where properties such as the 
Gaussian curvature also vary with position. 



APPENDIX B: GEOMETRIC POTENTIAL 

In this appendix we present two equivalent ways of de- 
termining the explicit form of the geometrical potential 
V{u.) valid for azimuthally symmetric surfaces like the 
Gaussian bump. The starting point is the general defini- 
tion introduced in Section ITTni 



F(u) = - fdA' G(u') r(u,u') , 



(Bl) 



with the Green's function F as defined in Ea. l|A10() . In 
the electrostatic analogy, V^(u) is thus the potential in- 
duced by a continuous distribution of "charge" repre- 
sented by the Gaussian curvature (with sign reversed). 
We shall derive an analogue of Gauss law for corrugated 
surfaces where the curvature of the surface (with sign 
reversed) plays the role of a continuous density of elec- 
trostatic charge. 

The first derivation makes use of the fact that V{u) is a 
scalar under conformal transformations. This symmetry 
can be checked explicitly by applying the same reasoning 
adopted for the equation satisfied by the Green's function 
in Appendix fXl to Ea. (jBl|) . In fact, upon operating on 
both sides of Ea. (|BH) with the covariant Laplacian and 
using Ea. ljAlOjI . the defining equation for the geometric 
potential can be cast into the differential form: 



D^DaViu) == G{u) 



(B2) 



The Gaussian curvature in Ea. (|B2l) can be written in 
conformal coordinates [31| as 

G{x,y) = -e-P^-^'VHdl + dlf-^ , (B3) 

where p(x, y) is the conformal factor introduced in Ap- 
pendix^ Similarly the left hand side of Ea. (|B2() can be 
expressed in conformal coordinates |3lj | as 

D-^DaVi^) = +e-P^^-yHdl + d^)Vix,y) . (B4) 



Upon substituting Equations IjBSp and HB4|) in Ea. l|B2|) . 
we conclude immediately that the geometric potential in 
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conformal coordinates V{x) is given by 

T./ ^ PJx, y) 

V{x,y) ^ — . 



(B5) 



Upon using Equations (|A6I) , ljA8|l and lA9p to substitute 
in Eq. (|B5|) , one obtains the explicit form of the geometric 
potential for the bump parameterized by the coordinates 



(r,0): 



Vir) 



f 

Jr 



dr 



,VW) 



(B6) 



where the a-dependent function l(r) was defined in 
Eq. (jA4|l and plotted in Fig. |2| The result of the integra- 
tion in Eq. (|B6|) is independent of (/) because of azimuthal 
symmetry. The upper limit of integration is chosen con- 
sistently with Ea. (jB2p as usually done in electrostatics. 
A second derivation of this result is obtained by mak- 
ing explicit use of the azimuthal symmetry of the bump 
and deriving a covariant form of Gauss' law which allows 
an intuitive understanding of the interaction between de- 
fects and curvature. This curved space version of Gauss' 
law illuminates the electrostatic analogy used throughout 
the text. The gradient of the geometric potential defines 
an "electric field" £"„ given by 



Ec, = -Dc,V{u) = fdA' G(u') D„r(u,u') 



(B7) 



where we have used Ea. (|Bl|) . One might expect that 
the flux of the vector E" through a closed loop is pro- 
portional to the enclosed Gaussian curvature in analogy 
with Gauss' law that relates the flux of the electric field 
to the electrostatic charge enclosed. To prove this asser- 
tion, we invoke the generalized Stokes' formula |2a| that 
relates the surface integral over A of the gradient of a 
field to its flux through the contour loop C: 



dA DaE" = - (t du°'-fjEp 



(B8) 

I A JC 

The covariant antisymmetric tensor 7^^ is given by 

lap = \^eQ/3 , (B9) 

where tap is the anti-symmetric tensor with e^^ = 
—e^r = 1- Similarly 7"^^ equals ^^ and the following 

identity holds 21] 



7 7<T/3 



COL 

^-6p . 



The tensor 7^'^ 



(BIO) 
laag'^ performs anti-clockwise rota- 



tions of ^ when acting on an arbitrary tangent vector V^, 
as can be checked by evaluating V^-^fa^Vp = jafjV'^V^ = 
0, where we have used the antisymmetry of jafs- Thus, 
the vector du^ja^ in Ea. (|B8|) represents an infinitesimal 
contour length times the inward unit vector perpendic- 
ular to it. The dot product with the field Ef^ then gen- 
erates the flux. To calculate the flux piercing a circular 



circuit centered on a Gaussian bump, we will need to 
explicitly evaluate 7^'", 



70 



(Bll) 



Upon using Ea. (|B7p we can rewrite the left hand side 
of Eq.(|B8l) as 



fdA Dc,E°' = fdA fdA' G(u') DpDf^Tiu, u') . 

(B12) 
If we now recall the defining Equation l|A10p of the 
Green's function r(u, u') and keep in mind that the 
Laplacian in Ea. l|B12|l operates on the variables labelled 
by u' (not u), we obtain using Eg. ljBSII a general result 
for the flux piercing a closed loop on the surface, namely 



du^jjE, 



d^u^ G(u) 



(B13) 



We can explicitely evaluate the right hand side of 
Ea. (|B13|l with the aid of Ea. ((T^ . By appealing to the 
cylindrical symmetry, in the special case of the Gaus- 
sian bump one can apply this covariant form of Gauss' 
theorem to find the radial field E^(r) in terms of the in- 
tegrated Gaussian curvature divided by the length of a 
boundary circle of radius r, with the result 



E'' 



1 



l{r) 



(B14) 



where we used Equations (|A2I) and (|B11|) . The angular 
component E'^ is zero everywhere. Note that E^{r) van- 

ishes linearly with r for small r and decays like re ''o 
for r ^ ro- From Ea. (|B14|) one can obtain the geometric 
potential V{r) by performing a line integral. 



V{t) 



dr' grr E^ 



dr- 



,,v^F)-i 



(B15) 



which matches the result previously obtained in Eq. IJB6|I . 



APPENDIX C: FREE ENERGY ON A 
CORRUGATED PLANE 

In this Appendix, we derive the effective free energy 
for a charge neutral configuration of defects confined on 
an infinite surface of varying Gaussian curvature with the 
topology of the plane. A general method was introduced 
in Ref.^(T| that allows treatments of the more compli- 
cated case of deformed spheres. A detailed treatment of 
boundary-effects is developed in Appendix ^ Here we 
simply assume that the size of the system is much larger 
than the size of the bump and that the boundary does not 
impose any topological constraint to the director of the 
liquid crystal. The results presented here match those 
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FIG. 17: (a) Defects with fixed core size a on a Gaussian 
bump encircled by a circular boundary of radius R denoted 
by B. (b) The size of the vortex cores varies with position 
when plotted in the conformal plane. One can avoid the sin- 
gularities associated with the defects cores by puncturing the 
conformal plane. This introduces circular boundaries d of 
varying radius at the position of each defect in the confor- 
mal plane, reflecting corresponding constant core radii on the 
Gaussian bump. 



obtained in Appendix IdI for free boundary conditions, as 
long as the defects are far from the boundary. Suppose 
that all the N^ defects have the same circular core radius 
a which does not depend on where they are located on 
the surface. This assumptions is justified if the radius 
of curvature — is much greater than a. In this limit, 
the microscopic physics that determines a is insensitive 
to the presence of the curvature and since the bump is 
locally fiat a is approximately constant everywhere. The 
starting point of our analysis is the free energy expressed 
in terms of the singular part of the bond angle ds (u) 



1 



Ka / dAg'^f'id^d, - A^){dpB, - Afi) . (CI) 



F 



The cores of the defects are excluded from the area in- 
tegral in Eq. (JC1|I , hence S' is a disconnected domain cor- 
responding to the corrugated surface punctured at the 
positions u^ = (?'i,'/'i) of the defects. In the conformal 
plane parameterized by the coordinate 5R(r) defined in 
Ea. (|A8|l . the defect cores are circles whose position de- 
pendent radius is given by aeX^"^^' . The boundaries of the 
"core-disks" are labelled by Q while the circular edge of 
the sample of radius R is denoted by B, see Fig. [T7I 
Upon introducing a "Cauchy conjugate" function x(u) 
defined by 

d^Bs -A^= -fjdpx , (C2) 

the free energy in Ea. llCl|) can be cast in the form: 

F^\kaJ dAg^f\d^xm3X) ■ (C3) 



In deriving Eg. ljCSp we used the identity 



which can be proved with the aid of Ea. (|B10|) and the 
discussion following it. Ea. (jC4p implies that the (covari- 
ant) dot product between two vectors after rotating each 



of them by ^ is equivalent to taking the dot product be- 
tween the two initial vectors. The integral in Ea. HC3() 
can be rewritten as 

Z_^lJdAD„{xD''x)-lldAxD^D^X, (C5) 



where 



Da.ixD'^x) = iyV9)MV9g'''xd0x) , (C6) 

and D^Dc is defined in Eg. ljAlU) . Upon taking an ad- 
ditional covariant derivative, we can recast Ea. ljC2p in 
the form of a Poisson equation for the electrostatic-like 
potential x(u)j 



D^D'^xi^) = -P(u) 



(C7) 



where the analogue of the electrostatic charge density 
p{u) is given by 



Y^ S{u,u,) 



(C8) 



It is useful to compare Ea. ljC7|l to Ea. (|B2p used to de- 
fine the geometric potential in Appendix B. Both expres- 
sions are Poisson equations, the only difference being that 
the source term of Ea. (|(y8p includes both the point-like 
charges of the defects and the Gaussian curvature with 
its sign reversed. Hence the Gauss law discussed in Ap- 
pendix ^ for the geometric field Ea applies also to i9qXj 
provided that Ea. (jB13|l is suitably modified to include 
the contribution from the topological charges of the de- 
fects: 



V9 



j du'^jjdpx = J d^n^ ( G(u) -Y^q, 

'^' (C9) 

where C is the contour enclosing an arbitrary surface A. 
This relation will be useful later. 

We can formally solve for x(u) in Ea. (|(]7p in terms of 
the Green's function r(u, u') found in Appendix 



X(u)= ( dA' p(u') r(u,u') 

J A 



(CIO) 



where boundary terms have been dropped using charge 
neutrality and the fact that the edge of the sample is 
assumed to be far away from the defects. The integral in 
Ea. (|C10|) can be evaluated, with the result 

X(u)=^g,;r(u,u,)- / dA' G(u') r(u,u') . (CU) 



(C4) Upon using Eg.HBljI. we obtain: 



Nd 



(u) = ^g,r(u,uO + y(u) 



(C12) 
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We note that the first term is singular at positions u^ but 
when X is substituted in Ea. (|C3|l the resulting energy is 
finite because the core of the defects are excluded from 
the domain of integration S. Upon substituting Eg. ljCTp 
in the second term of Eq. IjCSp we obtain: 

( dAxD^D^X - - f dAxp{u) 
Js Js 

= [ dAxG{u) , (C13) 
Js 

where we dropped terms involving the delta functions in 
Ea. (|C8|l because they vanish everywhere except at the 
coordinates of the defects which are excluded from the 
domain of integration S. Upon substituting Ea. HC12|l in 
Ea. (|C13|l we obtain: 

- / dAxD„D''x^y2<}^yi^^) 

+ f dA jdA' G(u) r(u,u') G(u') , (CM) 

where we used Ea. (|ljl() . 

To evaluate the first term in Eg. ljCSI) . we apply the 
generalized Stokes formula of Eo. ljBSp and convert the 
surface integrals into line integrals over the boundaries: 

„ Nd „ 

/ dA D^ixD^x) ^ Y.f du^XlJDpx 

du'^XloPx^SiX. (C15) 

B 

where the difference in sign between the two boundary 
integrals in Eg. ljClSj) is due to the fact that the outward 
normals for the paths Ci are oriented opposite to the 
normal for B, the outermost boundary of the system. 

To evaluate the last term in Ea. ljC15|) . we note that 
the flux through the distant boundary B due to a charge 
neutral distribution of defects is approximately zero, pro- 
vided that the integrated Gaussian curvature enclosed by 
the boundary is vanishingly small (see Eg. UCQjl ). Hence 



:(^du''-fjDpx^O, 



(C16) 



where we used the fact that x, defined in Eg. (|C2|l . is 
approximately constant on B since daO — A^ — 0. By 
contrast, the flux of drX piercing the boundary Ci in 
Eo. (jC15|l is approximately equal to the charge qi of the 
inclosed defect [i^- In evaluating the integrals around 
the inflnitesimal boundaries Ci, we used the fact that the 
function x(ui + a) evaluated on the "rim" of the defect 
core centered at u^ and of radius ae^^"') is dominated 
by a logarithmically diverging contribution due to the 
I*'' defect. This leading contribution is approximately 
constant on Ci. On the other hand, the non diverging 
part of x(ui + a) is multiplied by the perimeter of Ci and 



hence its contribution is of the order of a. The result of 
the integration will be insensitive to the orientation of 
the vectors along the boundary Ci, provided the defect 
core is small '46]. In this way, we find 

du'^loPxDaX - X(u» + a) * du'^j/drX 

~ q^x{u^ + a) . (C17) 

Upon substituting Equations HC17|I and (jC16p in 
Eg. (jC15|) . we obtain 

„ Nd 

/ dAD^ixD'^x) = Vg»x(u»+a) = 
Js 



Nd 



i=l 

Nd Nd 



2—1 -i— 1 j^i 



Nd 



+ ^(Z?r(ui,u, + a) , (C18) 



i=l 



where we used Ea. (jC12|) to substitute for x- Substitution 
of Eg. (1^18)1 and Eq.lEHJ in Eq.lEU then yields: 
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Nd 
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^EX1*^J^("»'"J 



(C19) 



= 1 j^i 



where we have used Ea. ljAlSj) to evaluate r(ui, u^ + a). 
The first term in Eg. ljC19|l is the free energy of the 
smooth defect- free texture (see Ea. H17|l and preceding 
discussion). The free energy difference -j^ between a 
charge neutral defect configuration and a smooth texture 
thus reads 

/MTf \ T Nd Nd 77 ^d 

^ 9Z^Z^'?»9ira(r,,(/)„rj,0j)-|--— 2^gi 



Ka 



«=1 i#i 



Ka 



Nd 



(C20) 



where the subscript a in the Green's function indicates 

-2^^cos(0-0')].(C21) 
a a 

In order to absorb the core radius a in the inter-defect 
interaction, we used the elementary algebraic identity 



Nd Nd Nd 

1=1 i—1 j^i 



(C22) 



valid for charge neutral configurations. The core energy 
Ec in Eg. ljC20|) was added by hand and represents short 
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FIG. 18: (a) Schematic illustration of the boundary direc- 
tor texture corresponding to free boundary conditions. The 
vector order parameter orientation close to the edge of the 
sample does not vary appreciably as one moves along the ra- 
dial direction and it is parallel to itself at every point on the 
boundary (b) Tangential boundary conditions. The vector or- 
der parameter is locally aligned to a wall located at the edge 
of the sample. 



distance physics on scales less than or equal to the core 
radius. Although the second part of the last term in 
Ea. (|C20|l arises as a position dependent self energy, it 
has the same functional form as the geometrical poten- 
tial discussed in Appendix IbI and hence depends on the 
global shape of the surface. 



APPENDIX D: BUMP WITH A BOUNDARY 

The aim of this appendix is to study the energetics of a 
singular vector field on a bumpy surface of circular shape 
and finite size R. To evaluate the ground state energy 
in Ec ljCip . we first need to solve the covariant Laplace 
equation for the bond angle field 9{u) in the presence 
of Nd defects at positions u,. This is more easily done 
by switching from 6{u) to the conjugate field x(u), as 
shown in Eq. HC3|) , and solving the Poisson Equation (jC7p 
satisfied by x for both free and fixed boundary conditions 
(see Fig. tB- 

The bond angle field satisfies free boundary conditions 
if the following relation holds on the circular edge B: 

a,0|,=fl, - , (Dl) 

while for fixed tangential boundary conditions we have: 



dd,0\r=R — 



(D2) 



To understand Ea. HU2|l . recall that we measure the bond 
angle 9 with respect to a rotating basis vector Er in the 
radial direction. With this convention, 6 is equal to a 
constant when the vector order parameter is aligned with 
the circular boundary B. We can convert Equations ljDl|) 
and ljD2|) into boundary conditions to be satisfied by the 
conjugate field x O'n. B. Upon substituting Ea. (jDl|l in 
Ea. (jC2|l and using the fact that Ar is equal to zero, we 
obtain the constraint that x(u) satisfies on B in the case 
of free boundary conditions: 



This corresponds to a Dirichlet problem where x^(u) 
evaluated on the boundary B assumes an arbitrary con- 
stant value, c: 



X"(S) 



(D4) 



Upon substituting Ea. HD2|) in Ea. (|C2p . we obtain 



drX 



A^ 
70'' 



(D5) 



since ^t ~ 0. Upon substituting Equations HBlip and 
© in Ea. (|D5|l we obtain the boundary condition on x 
that corresponds to Ea. ljD2|l 



drX 



N\ 



\r=R 



1 

R ' 



(D6) 



where the superscript indicates that this is a Neumann 
boundary problem with the normal derivative assuming 
a constant value. 

To solve the Poisson Eq. ljC7p with Neumann or Dirich- 
let's boundary conditions in terms of suitable Green's 
functions we exploit a covariant version of Green's the- 
orem expressed in terms of two invariant functions of 
position '(/'(u) and (p{u) 38]: 

dA {(p{u)DaD'^ij{u)-^{u)DaD°'Lp{u))^ 
Js 

du"7„'3(^(u)a^V(u)-^(u)a^(^(u)) . (D7) 

By applying Ea. (|D7p to f{u) — x(u) and ?/'(u) — 
r(u,u') and using Equations (|A10|) and l|G7|l we obtain 

X(u) = fdA' r(u',u) p(u') 
Js 

- jdu"^^jT{x^,n)d'px{^), (D8) 

where u and u' have been exchanged [43 . The boundary 
conditions for the Green's function can be conveniently 
chosen to eliminate unknown quantities in Eq. IJD7|) as in 
fiat space |3a |. 

For the Dirichlet 's problem, we choose the Green's 
function F^ so that it vanishes when u' is on the bound- 
ary B 



F^(B,u) = 



(D9) 



Upon substituting Eg. ljDQp in Eq.jDS} and noting that 
x(u') is constant on the boundary B (see Ea. (ID4ll ). we 
obtain: 



X^(u) = (dJ^ F^(u',u) p(u') 
Js 



d4,x\r=R = 



(D3) 



+ x'^iB) i du'"7/5^r^(u',u) , (DIO) 
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The contour integral in the second term of Eq. (|D10|) cor- 
responds to the flux piercing B which is, in turn, equal 
to the (unit) charge of the singularity located at u |48j . 
The final result reads 

X^(u)= jdA' p(u') r^(u',u)+x^(S) , (Dll) 
Js 

where x^(i?) can be set to zero, since the energy in 
Ea. ()C3|l is only defined in terms of derivatives of x- O^ie 
can check that X'°(u) in Ea. (|Dll|) satisfies both the Pois- 
son's Equation l|C7|l and the required boundary condition 
in Ea. HD3|l . This can be more easily proved by noting 
that r-°(u',u) is symmetric under exchange of its argu- 
ments u' and u p9l |. 

A similar reasoning applies to x^(u)- However, we 
cannot choose the Green's function so that the second 
term of Ea. (|D8|l (that contains the unknown quantity 
x(u')) vanishes. In fact, by invoking Stokes theorem (see 
Ea. (|B8|) ). we note that 

(f du'°'-fJd'fjT{u',u) = l, (D12) 

where ja {'^' = B) is constant on the circular boundary 
B and can be brought out of the integral. An appropriate 
choice of boundary condition on T^ that satisfies the 
constraint in Ea. l|D12|l is 



a,,r^^(r',0';r,'^)|.' 



1 



27r7„/3(i?) 



2nR 



, (D13) 



where we used Ea. (|Bll|) . The a-dependent function l{r') 
was defined in Ea. l|A4|l . Note that 1{R) c± 1, f or i? > ro 
(see Fig. |2). By substituting Equations HD13|) and HD6|) 
in Ea. ljDSp we obtain 



X^(u) = fdA' r^(u',u) p(u') 
Js 



2^./o 



d0'x(i?,0') 



VW) 



2Tr 



i'tN 



r^(i?,0';r,</)), (D14) 



The last two integral are constant and hence can be 
dropped. We can check explicitly that x^(u) satisfies 
Ea. (jD6p by evaluating the radial derivative of x^(u) in 
Ea. (|DT4JI 

drx''{r,^)\r=R= fdA' p(u') 9,r^(i?,<^';r,0)|.=i?, 
Js 

(D15) 
where the radial derivative of the Green's function as- 
sumes the constant value derived in Ea. (|D13|) . provided 



that r^(u',u) is constructed so that it is symmetric un- 
der exchange of its arguments u' and u (see Ea. (|D19|) '). 
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FIG. 19: Schematic illustration of the method of images. The 
image defect is of the same sign for free boundary conditions 
(a) and opposite for fixed boundary conditions (b). Defects 
closer to the center of the circle have images further away 
from it. 



With the aid of Equations (|G8p and ifT^ , we obtain 

Upon substituting Equations ljD16|l and (|D13|I in 
Ea. (jD15|) . we conclude that the the Neumann boundary 
condition in Ea. ljD6|) is satisfied provided that 



^ (J, = 27r . 



(D17) 



It is reassuring that the topological constraint on the 
vorticity of the field imposed by the presence of the wall 
arises as a natural requirement within this formalism. 
Similarly, the Poisson's equation for x^(u) is automati- 
cally satisfied. 

We are now left with the task of guessing the Green's 
functions for the Dirichlet's and Neumann problems sat- 
isfying the boundary conditions in Equations l|D9|) and 
(|D13|I respectively. In both cases the Green's function 
can be determined by the method of images applied in 
the conformal plane. For every defect with radial coor- 
dinate ri we need an image defect of opposite (equal) 
topological charge at position r^ to ensure that Dirich- 
let's (Neumann) boundary conditions are enforced (see 
Fig. I19|l . The radial coordinate of the image defect r^ is 
determined by the relation: 



Wi) 



5R(r,) 



(D18) 



Except for the coordinate change r — + 5ft(r'), a similar 
relation arises in elementary electrostatic problems in 
flat space Q. A geometric argument that justifies this 
choice of images in flat space is illustrated in Fig. 1201 

Once the position of the source is chosen according to 
Eq. (|D18|) , we can express the two Green's functions with 
the concise notation r^/^ as follows: 
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r^/^(u,u') = --^ln[SR(r)2+3fi(r')'-23fi(r)3?(r')cos( 

47r 



]±-ln[5i(r)^- 






-2di{r 



njRY 



cos 



(D19) 



where we have introduced a function /(?'') to make 
r-D/ATj-^^ u') symmetric under exchange of its arguments 
and to remove a singularity at r' = 0. Note that we 
can add f{r') since the defining equation of the Green's 
function docs not contain derivatives of r', only of r. 



/(.') = ^ln 



^{r') 



^{R) 



(D20) 



The plus and minus signs in Ea. (|D19|l insure that the 
Dirichlet and Neumann boundary conditions respectively 
are obeyed. In what follows the sign placed above in the 
symbols ± or =F always indicates the choice suitable for 
the Dirichlet 's problem while the one below refers to Neu- 
mann's boundary conditions. One can explicitly check 
by substitution that the symmetrized Green's functions 
r^/^(u, u') satisfy the correct boundary conditions, as 
long as the plus sign is chosen when F^ is substituted in 
Ea. (jD9p and the minus sign when T^ is substituted in 
Ea. (|D13|l . Note that, without the extra term /(r') in the 
expressions for both Green's functions, F^ would not be 
equal to zero on the boundary B and the last term in 
Ea. (|D14|) would not be constant when T^ is substituted 
in. 

Once the Green's function is obtained, one can readily 
write down x^ (u) by dropping the constant terms in 
Equations IjDlip and HD14|I 



D/N 



Na 



(u) = J2q,T^/^ {u,n,) 



- fdA' G(u') r^/^(u,u') . (D21) 
Js 

The Gaussian curvature is given by the covariant Lapla- 
cian of the geometric potential introduced in Ea. ljB2p . 
Upon integrating by parts twice the second term in 
Ea. (|D21|) and applying Stokes theorem repeatedly, we 
find 



Na 



D/N 



(u)-^g,r^/^(u,uO + y(u), (D22) 



where we assume that i? 3> tq so that we can neglect 
boundary terms. The geometric potential V^(u) has the 
same functional form previously discussed in Appendix 
IbI despite the change in the Green's function. 

The evaluation of the energy stored in the field now 
proceeds along the lines sketched in Appendix El with 
the only caveat that one needs to choose the appropriate 
Green's function in Ea. (|D19|l . In the case of Dirichlet's 
boundary conditions one can prove that the boundary in- 
tegral in Eq. ljG16p still vanishes by virtue of the fact that 




p, X 



FIG. 20: A topological defect located at position P in a circu- 
lar domain of radius OQ in flat space. Fixed boundary con- 
ditions are obtained by placing an image defect of the same 
sign at a distance OP' from the center such that OP OP' = 
OQ . The two triangles AOQP and AOQP' are similar and 
ZOQP = ZOP'Q = TV - e'. By the theorem of the external 
angle, we conclude that 6' + 6 = (j) + n as long as Q lies on the 
circumference B. This is equivalent to the boundary condition 
in Ea. lD21 if a non rotating vector basis is used. Similarly, 
for free boundary conditions the symmetric Green' function 
evaluated is constant on the boundary if the image defect is 
negative. Since =^ = 2£, the potential In(PQ) - In(P'Q) 

(generated by the defect at distance OP and its image) is 
constant on the circumference of radius OQ. 



X is constant on the boundary B and the defects configu- 
ration is charge neutral. For the Neumann's problem we 
have: 



,N 



j> du^jjDpx^ 



47rln(3ff(i?)) 



(D23) 



where we assumed that R ^ rg. In this limit ^{R) is 
approximately equal to R, as can be checked with the 
aid of Equations (|A8(I and (|A9(I . 

All the remaining intermediate steps to derive the free 
energy follow as in Appendix [CI without further assump- 
tions. We can readily generalize Ea. l|G20|) to evaluate 
the energy stored in the singular field in the presence of 
a boundary in the case of both free and fixed boundary 
conditions. We assume R ^ ro, but the defects do not 
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need to be far away from the boundary. The result is 
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(D24) 



where Fq is defined in Ea. p7|) . The Green's function 
expressed in scaled coordinates reads 



i yX{ , Xj j 



47r 



In [xf + x^j — 2xiXj cos (0i — (j)j)] 



±— In [xfx'^ + 1 — 2xiXj cos {(pi — 0j)] . 

(D25) 

In the case of Neumann's boundary conditions, we have 
suppressed a term diverging like ln(9fi(i?)) associated with 
the boundary contribution in Eg . (10231) . Ea. l|D25|) is ex- 
pressed in terms of a dimensionless defect position Xi 






(D26) 



The plus sign in Equations ljD25p and (|D24p is to be cho- 
sen for Dirichlet's boundary conditions and the minus 
sign for Neumann's. The last term in Ea. ljD24|l repre- 



sents the interaction 

located at Xi and the boundary 



C/jj ' (xi), between a single defect 



U, 



D/N 



(x,) = ±K, 



-E' 

i=l 



in 



In [1-^2] . (D27) 



Note that the q-dependent prefactors of U^ (xi) and 
the quadratic correction to the curvature interaction (III 
term in Ea. (|D24|l have the same magnitude. This is not 
a coincidence but a clue to their common origin. As the 
geometry of a plane is modified, either by creating a vary- 
ing curvature or imposing boundaries, the defects feel an 
additional interaction caused by the conformal transfor- 
mation of the underlying space. This line of reasoning is 
powerful and it has been pursued in Ref. 30l| to explain 
some basic features of the interaction between defects 
and curvature without explicit recourse to the Green's 
function techniques adopted in this work. 
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The integrated Gaussian curvature in the microscopic 
disk is vanishingly small. 

If the defect core is very large, it may be necessary to 
place images within the defect core itself to impose a 
desired boundary condition on its rim. This is not the 
regime considered in the present work (see Ref. |4H| for 
similar calculations performed in fiat space). 
These manipulations are common in electromagnetism, 
see for example Ref. |39|. 

This assertion can be proved by applying Stokes Theorem 
(as stated in Ea. llB8ll with Ep replaced by c)^r(u',u)) 
and using Eg. lAlOt 'l to evaluate the surface integral. 
This assertion can be proved by applying Green's theo- 
rem in EqinHto V(u) = r(u,u') and ipivi) = r(u, u") 
and noticing that the right hand side vanishes if the 
boundary condition in Ea. llD9t is assumed. We can then 
conclude that r°(u',u") = V°{\i",\i') in analogy with 
familiar results in 3D electrostatics 39). 



